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Abstract 

In the framework of the Kompaneets approximation the propaga- 
tion of a shock front (SF) in inhomogeneous medium with the power- 
law density decrease at the exponent n = 2 is investigated. For this 
important case corresponding to outer regions of the solar and stellar 
coronas an unexpectedly simple exact solution clearing a structure of 
the general solution for an arbitrary monotonic density of medium 
is found. Our results for plane-layered medium are compared to the 
Korycansky one for off-centre explosion in radially stratified medium. 
The relation between the solutions in these media are found and a 
new exact solution for a noncentral explosion in the case of density 
singularity on the finite radius is received. 
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Introduction 



As really seen from the well-known Sedov-Tailor formula in the homogeneous 
medium [1, 2] 

R = [E t 2 1 pop (1) 

for the dependence of the radius R of a shock front on time t from the 
moment of explosion, the explosion energy E and the medium density p , the 
SF would move faster towards the direction of density decrease if explosion 
occurs in inhomogeneous medium. However, for inhomogeneous medium 
the problem becomes more complicated, being not self-similar, and requires 
significantly more complicated calculations even for computational analysis 
(see the review and papers [3, 4, 5] and the bibliography therein). 

Assuming the uniform pressure along the envelope behind the SF 
A. S. Kompaneets proposed the equation [6] (see Eq. (3) below) for strong 
explosion in an exponential atmosphere. In that case he was showed that 
the SF was accelerated and could get infinity for a finite time. (Surely, the 
approximation fails earlier and account for a non zero density at infinity 
[7] restricts values of the SF velocity too, i.e., there is no controversy with 
physics or common sense.) For densities decreasing not so fast, namely, for 
the power-law density dependence p(z) oc z~ n , often met in astrophysical 
objects, the behaviour of the SF is defined by the exponent n [2, 8, 9]. The 
asymptotically exact solution for the head part of the SF moving towards the 
direction of density decrease indicates that the Sedov-Tailor SF deceleration 
is down when the density gradient (or n) is up 1 . Particularly, it is confirmed 
by the exact solution of the Kompaneets equation discussed below for n = 2 
[12] corresponding to the wind regions in outer parts of the solar and stel- 
lar coronas (however, with the wind motion not taken into account). This 
solution having an extremely simple form is analysed in this publication. 

Following Kompaneets [6], instead of the real time t we define the new 
temporal variable y as 

dy = dt[E \( 1 2 -l)/(2p V(t))] 1/2 , (2) 

including into it the inverse square root of the time dependent volume V(t) = 
7r / r 2 dz. After substitution the coefficients of the equation for the SF r = 

1 Note that at n = 3 the SF moves at a constant velocity. For n greater than 3 the 
SF is accelerated and since n = 5 it gets infinity for a finite time [8, 9, 11]. The solution 
qualitatively becomes similar to the Kompaneets one in an exponential atmosphere. 
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r(y, z) do not depend on time any more but only on the coordinate z. Along 
this coordinate the density changes as p{z) = po/(z/z ) 2 . Here po is the 
unperturbed density in the point of explosion z . Some another relations 
important for applications but not for the current formal consideration are 
presented in Appendix 1. The Kompaneets equation for the shock front can 
be written as: 

(dr/dy) 2 <p(z) - [{dr/dzf + 1] = 0. (3) 

Here <p(z) = p(z)/p . In the case of n = 2 the solution could be immediately 
written, however, we will consider it as a particular case of the general sit- 
uation to clarify its structure in more general cases of astrophysical interest 
(change of the exponent n from 16 to 6 and 2 in the solar corona; n = 3 in 
active nuclei of galaxies). 

The plan of this paper is the following. In the beginning we give the 
quadrature parametric solution of the Kompaneets equation (3) (see [6]) for 
an arbitrary density distribution of the medium using a supplementary func- 
tion £(y,z) and an arbitrary function //(£). Regarding the medium at the 
initial stage of expansion as homogeneous we conclude from comparison to 
the Sedov solution that p is to be zero (see [6]). At p — for the inverse- 
square law of the density decrease we find an explicit solution in two regions 
adjoined to the opposite leading points of the SF. In this two regions the 
solution consists of two non-adjoined segments of the same sphere. In the 
intermediate region the function p is not zero as well as in the Silich & 
Fomin case [7], where the density decreases exponentially to some constant 
value. The solution obtained in the intermediate region complements the 
SF to a full sphere. It is shown that at short times (y — > 0) the interme- 
diate region disappears. The transformation is also found of this solution 
in a plane-layered medium to that for an off-centre explosion in the radi- 
ally stratified medium with the density distribution having a singularity at 
a sphere surface with a finite radius, i.e., one more new solution is obtained. 
(Preliminary comparison is carried out to the Korycansky solution [8] for 
a power-law radially-stratified medium). An extremely simple form of the 
solution obtained admits explicit calculations which can be important for 
general analysis of propagation of strongly non-linear waves in an inhomo- 
geneous media. A square-root singularity at the leading points is revealed 
as an inherent feature of £ for an arbitrary density dependence on coordi- 
nates. This feature allows some conclusions to be made about SF properties. 
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The fact is that different representations of the solution is required in dif- 
ferent regions. Probably, in the vicinity of every leading point we deal with 
the only non-linear wave whereas contribution of both of these waves (non- 
linear "interference") should be matched by choosing the function \i in the 
intermediate region. 



Initial condition 

The complete integral of Eq. (3) dependent on two arbitrary constants (ac- 
cording to the number of independent variables) can be found via separation 
of variables (see, for example, [13]). Regarding the integration constant \i as 
an arbitrary function of the separation constant £ and substituting the latter 
with the function of two independent variables z), which can be derived 
from the condition [13] 

dr/d£ = 0, (4) 
we find the following expression for the general integral: 

r = - sign(^ - zq) I dxJ^ip(x) - 1 + £y + //(£)• (5) 

Jz 

This general integral and the condition (3) for the function £ give its implicit 
expression via z and y 

y = sign(z - zq) / dx ==== = - fj, , — = // , (6) 

zo yi <p(x) - 1 ^ 

from where for the solution describing the SF we have 

rz d x 

r = sign(z - z ) / —====== +fi- (7) 



At y — > 0, z — > Zq, r — > 0, 



1 i 

===\z-Zq\ + //-£//, yK-===\z-Zo\-fi', (8) 



or 

V ^-^o / f 2 -l' \Z-Zq / 



(9) 
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Generally, the supplementary functions [i and £ are considered not expanded. 
Excluding £ 2 gives: 

It is obvious that approaching the Sedov limit 

r 2 + (z - z ) 2 = y 2 (11) 

is possible only if (disregarding possible though non-essential shift of initial 
moment of explosion) 

^^-0, ^-0, y^O. (12) 

r y 

Finally, for this limiting case of medium homogeneous at small distances the 
following relations occur 

fi' = 0, fl = 0. 

We will see below that the situation is more complicated, however, for while 
we can use the relations derived for the exponent of interest n = 2. 

Solution for n = 2 in the vicinity of leading points 

As mentioned above, this case is important for the astrophysical applications 
(outer parts of coronas, the Solar wind region, etc.). However, it also gives us 
a clear example of how to express the solution via supplementary functions. 
At ip(z) = 1/z 2 (zq = 1) the integrals are expressed via elementary functions 
but, as seen, separate regions in the z,y plane should be considered to pro- 
vide positiveness of the expressions under the square root. Let z satisfy the 
inequality z 2 > z > z cx . In this region adjacent to the leading point of the 
SF we obtain (fj, — 0): 

fV > 1, f Lx= z cx , (z = 1), 



e-i-^je-z 2 , (*>i), (is) 

z 2 (l-e 2 y) 2 



Ae 2 y(ey - z)(z- e~v) 



(14) 
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Figure 1: The sphere SF corresponding to strong explosion in a plane-layered 
medium with the inverse square density decrease, a) The half-sphere adjoined 
to the upper leading point 22 (t) moving towards the density decrease, b) The 
surface of the sphere segment adjoined to the lower leading point Z\{t). c) The 
surface of the spherical layer complementing the SF to the sphere. This part 
of the solution contrary to a) and b) corresponds to fi 7^ (the region of the 
"interference" of waves). The parts of the SF corresponding to the solutions 
in different regions (see text) are shifted here each off another. 



Here 

z 2 = e y , Zl =e~ y , z ex = coshy, r \ Zcx = sinhy. (15) 

Calculating 



[,(1 + ^-2^ 

q 4e 2 y{ e y -z){z-e-y) ' 1 ; 



and 

P 2 _ 2 - 

" 4e 2 y(ey - z)(z- e _ ») 



e _ z2= z\l + e 2 y-2zeyf 
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we find that 

Z>Z CX , ^l-Z 2 = , Z{Z ~" CX) >0, eU x =^ex- (18) 

yJ(Z2-z)(z-Z{) 

V^ = T( ZZcx ~ 1 f . ei^=i- (19) 

Substituting these relations into Eq. (13) we get an extremely simple solution: 

r = ^/(^ - z) (z - zi) , z 2 > z> z ex . (20) 

Apparently, it is the equation of the hemisphere hanging on the upper leading 
point z 2 with the centre at z ex and radius (z 2 — z\)/2 (Fig. la). Note that in 
this region really £ > z (17). At z 1 < z < z~ x < 1 with £ 2 from (14) we have 



= y/e~z i -yfe Z l, (21) 



^ 2 -z*= Z{Zcx Z) , y / e T T= , 1 ZZcx , (22) 

(z 2 -z)(z-z l ) y (z 2 - z)(z - Zx) 

and, finally, (Fig. lb) the solution reduced to the same form 



(z 2 - z)(z- z^, zz cx < 1 

describes the segment surface of the same sphere leaning against the lower 
leading point z±. The square root singularities of £ at the leading points are 
the general property for an arbitrary n substantially used in [9]. 



Intermediate region 

Here we return to the initial general form of the solution with the arbitrary 
function which is found from matching conditions at the boundary between 
the regions (compare to Silich & Fomin [7]) 

r Lx-0= r Lx+O, (^ex = £ Lx) (23) 

MO = 2 V / e T l - 2£ In {£ + Je^l) (24) 
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Figure 2: The supplementary functions £ and /i parameterising the solution 
for the SF. The square-root singularities of £ at the leading points arise for 
an arbitrary n. 
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/ i'(o = -2in{e+ v /e r T} 

Substituting these relations into 







dr 



i; 



£(p dx 



we get the equation for £ at z cx > z > z~ x 



+ y + /AO 



y = In 



Using the definition 
we obtain 



£ = cosh a 



a = In , 



z 2 (zz 2 - 1) 



\ (*2 -Z) 

« l*«= « L- X i= °- 

Thus, the value £ (Fig. 2a) in this region is given by the relation 

£ Ux= coshy = z cx , £ | -i= 1. 



£ = cosh In 



ev - z 



i y~ z y- z cx 



It is seen that 



ex ^ex 



|"^0, jj! ^ (z = Z cx ). 



(25) 



(26) 



(27) 

(28) 
(29) 



(30) 
(31) 

(32) 



Therefore /x is continuous at the lower boundary but together with the first 
derivative it jumps at the upper one (Fig. 2b). 

M&x) = 2(sinhy - t/coshy), //(£ cx ) = -2y; £ cx = £(z ex ). (33) 

Hence, in the region z cx > z > z~^ the expression can be written for £ 

I I \ev{e*z - 1) c» -z 



e» - z 



+ 



ey(eyz -i)f' 



(34) 



or 

zsinhy 

( = ^ -*)(,-„) • (35) 



The following relations also can be found 

z(z 2 - zi) 



2^/(^2 - z)(i - ^) 



(36) 



^ 2 _ 1 = (^coshy- l) 2 ^ 2 _ ^ = z 2 (cosliy - zf 
(z 2 -z)(z-z i y ' (z 2 - z)(z - Zi)' 

Finally, in the intermediate region (Fig. lc) the following formula is valid 



r = ^ e _ Z 2 + ^rr l (38) 

and according to (37), 

r = ^(z 2 - z )( z - Zl ). (39) 

Thus, we get the complementary surface of a spherical layer, which together 
with another segments above forms a complete sphere. It is a complete and 
extremely simple solution of the problem. We also present it in an explicit 
form indicating temporal dependencies of the envelope radius and centre: 

r = z [(exp{y/z } - z/z )(z/z - exp{-y/ z })} 1/2 . (40) 

The solution can be also presented in the canonical form: 

r 2 + (z - z cosh(y/z )) 2 = z 2 sinh 2 (y/,2o)- (41) 

At every moment the solution is a sphere with the centre at z = z cosh(y/,2o) 
and the radius growing as z sinh(y/zo). The dependence of the real time t 
on the Kompaneets variable y 



t oc J dy sinh 3 / 2 y 



can not be exactly expressed via elementary functions even for this simplest 
case but the asymptotical relation is obviously simpler. Such a perfect form of 
the solution in an inhomogeneous medium is very surprising. It is much closer 
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to the Sedov solution than to the Kompaneets one in full accordance with 
the comparatively slow density dependence on the coordinate (see footnote 
1). Nevertheless, the velocity decrease of the leading point at large time 
(z 2 ^> Zq)) dz2/dt oc t -1 / 3 is slower compared to the homogeneous medium 
dz 2 /dt oc t~ 3/5 . Note that at large y > 1 



z cx ~ 2 ' Z ex 



2z 1 . 



At small y < 1 



-i 



2 



When approaching to the point of explosion, the intermediate region gets 
narrow and disappears at all (Fig. 3). 



Comparison to results of the conform transformation 
for radially stratified medium 

Here we compare the results discussed above to the brilliant solution by 
Korycansky [8] for off-centre explosion in a radially stratified medium. In 
the Kompaneets approximation the equation x(i?, y) = 9 for the SF can be 
written as: 




1 




+ 



R 2 



^{R) = 



p(R) 

Po 



(42) 



Here R, 9 are the spherical coordinates, R — a, 9 — 0, (y — 0) are the coor- 



dinates of the explosion. At p(R) 

(f> = ax, a 



Pol 



\ the substitution [8] ( 



leads to the equation in an effective homogeneous medium 




1 

Vc 




+ 



2a 



n\ 



(43) 



The solution of [8] is given by the relations: 

2Ccos0= l + ( 2 -x 2 , 



_ y 

x = —, 

Vc 



(44) 
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Figure 3: The dependence of specific points of the SF on "time" y. It is seen 
that the region /i 7^ vanishes at y — *■ 0. At every moment the solution is a 
sphere with the centre at z — cosh y and radius growing as sinh y (zq — 1 is 
accepted). 



2-n 



cos ■ 



(2-n)0 



(IT- 


~y(2-n)~ 




2a 



(45) 



For n = 2 the following special solution exists [8]: 

a Va/ 

It can be easily found using the transformation from the spherical case 



(46) 



9y, 



R 2 



dR 



(47) 
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to the effective plane one with the help of the substitution: r] = ln_R, y — > ya: 



2 / o \ 2 



The same result can be obtained from the general formula (45) in the limit: 

Ra ^-n)/2 ^ exp t^l i n * _> 
2 a 

, 2-nR 1 /2-n\\ 2 R . . K } 

It can be easily verified that main parameters of the surface (46) are 
very close to that of the sphere discussed above. Correspondence to the 
case of homogeneous medium can be verified using the transformation to the 
cylindrical coordinates 



pY 2 

1 =tan 2 



,z 

At small time 

^2 („.U\2 l„2 



V_y _ ln 2 V/^jK*^ p = R s [ n e, z = Rcos9. (50) 

z z 



^ W*)-**'* - p (51) 



,zj 1 + \nz/ ' zq 

it becomes the Sedov solution in ordinary variables: 

p 2 + ( Z - Zo y = y 2 . (52) 

Correspondence of solutions for plane-layered and radi- 
ally stratified media and a new solution for SF in the 
case of density singularity on the finite radius 

The Korycansky transformations can be generalised. The following substi- 
tution of variables 

— = ln— (_R = aexp — ) (53) 

Zq a Zq 

transforms the Kompaneets equation for r(z,y) in a plane-layered medium 
(3) (for cylindrical coordinates z, r) to that for y) in a radially stratified 
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medium (42) (for spherical coordinates R,x)- This substitution also trans- 
forms the density dependence tp(z) in a plane-layered medium to ip(R) in 
spherical one according to the relation 

HR)=^o\n^ 2 (54) 

and vice versa: 

R 2 z 1z 

(p(z) = tp(R) — = (a/z ) 2, ip(aexp — ) exp — . (55) 

Z Zq Zq 

Corresponding solutions of the Kompaneets equations are related as 

r = z x- (56) 
Note that in this case the exponential density dependencies 

<p(z) = exp{-/?z} (57) 
are transformed to the power ones 

il>(R) = ^ exp{-[3z In -} oc R-W«>+ 2 ) . (58) 

The equivalent height of the atmosphere for the plane case determines the 
exponent for the power law in the radially-stratified medium. Particularly, 
this fact clarifies existence of the solution sequence in the Korycansky case 
[8], which therefore are connected with Kompaneets solutions [6]. 
In the special case 

(p(z) oc 1/z 2 

the exact solution for which is constructed above (see (40)) the law 

ip(R) oc l/R 2 \n 2 (R/a) (59) 

gives the density dependence very close to the inverse-square law at large 
distances and having a singularity at the surface of the sphere with the 
radius a. 

Correspondingly, there is a new formal solution with a fixed point of 
explosion R = ae: 

X(R) = ^{HR/a) - Zl /z )(z 2 /z - \n(R/a)), (60) 

where 21,2 (y)/£o are defined by the formula (15). 

The substitution (53) is also suitable for computations relating to the 
off-centre explosions in radially stratified media. 
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Conclusion 

The results obtained are in good agreement with the known self-similar so- 
lutions [14, 15, 16]. These results can be useful not only for the Solar corona 
but for another astrophysical situations (compare to [17, 18]) where the in- 
homogeneity of media is important and, moreover, the inverse-square law 
of the density decrease occurs very often 2 , for example, for Supernovae in 
the periphery of gigantic clouds, for outbursts in nebulae [20], for varieties 
of galactic fountains [21], for galaxies with star formation outbreaks [3] and 
also for more exotic cases such as the Supernova in the NGC 6946 nearby 
another Supernova (i.e., in its wind region) [22] recently discovered by the 
Hubble Space Telescope or for galaxy collisions. 

Particularly, the off-centre location of pulsars compared to Supernova 
Remnants can be naturally explained by displacement of the envelope to- 
wards the direction of the density decrease without attracting the idea of the 
non-symmetrical explosion and recoil. The solution in the medium with the 
singularity of the density at a sphere surface of a finite radius describes the 
bypass of the photosphere by the SF and is of great interest for the Sun and 
stars. The asymptotical solutions for the leading point reveal possible non- 
monotonic behaviour of the SF velocity at the monotonic density decrease 
(see [10]). The solutions obtained admits an interesting generalisation if the 
constant energy E is replaced by the function of time E(t) which allows to 
take into account losses [21] as well as an energy supply (see notice in [10]). 
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Appendix 1 



The Kompaneets equation for the SF is derived by equating the nor- 
mal component of the shock front velocity D = (p/(p — p 2 / p')) 1 ^ 2 ob- 
tained from boundary conditions to its formal one —(df/dt)/\Wf\. Here 
p = (7 — 1)XE /V is the pressure behind the shock front, 7 is the adiabatic 
exponent, E is the energy of the explosion, p and p' are densities ahead and 
behind the shock front, p'/p = (7 + l)/(7 — 1), V(t) is the volume of the 
envelope V(t) = n J r 2 dz, with the integral taken over the part of z axis, 
restricted by the envelope. The constant A, is of about to 2-3 and accounts 
the deviation of the pressure near the shock front from the average value, 
/(r, z, t) =0 is the shock front equation. 



Appendix 2 

Some notes should be made about the leading point behaviour for ip(z) = 
(z/z )~ n at n differing from 2. Relations (5), (6) and conditions r(zi t2 ) = 
of intersection of axis z by the SF give the following equations for the values 

y = jT dz'(<p(z')) 1/2 = 11 dz'( V (z')) 1/2 . (61) 

In the vicinity of the point z 2 the value £ is large and an expansion by l/£ 1 
results in the following expression for £ 2 [9]: 

e » [2 (z 2 - z) (y^ 2 )) 1/2 l 1 l Z2 du ( V (u))- 1/2 (62) 

Correspondingly, 

r » {l/£(z, y)) ■ fdu (<p(u)r 1/2 , z - z 2 (y). (63) 

Thus, the function £ has a square-root singularity with accordance to (14). 
These relations can be used to estimate the volume V(t) required to describe 
motion of the SF in terms of real time t. For z 2 we obtain: 

Z2it) = h • ( >>2/(n-5) ' ( 64 ) 

(1 - sign(n - h)t/U) n ' 



U = Jl6irzlp /(n + 2)(n- 5) 2 E A(7 2 - 1). (65) 



16 



As it follows from the relation (61) in the case of n > 5 [8, 9] "disastrous" 
acceleration takes place for the mentioned part of the SF which approaches 
the infinity z 2 = oo for a finite time. The expression for the velocity can be 
written as 

dz 2 /dtoc (U -t)(™- 3 )/( 5 -™). (66) 
The specific time t* of "disastrous" acceleration is 

U » (z 5 oPo /E ) 1/2 (67) 

(see also [9]). Despite the fact that we used here the solution valid only in the 
vicinity of the leading point the relations (64) and (66) are enough adequate 
to the phenomena. However, the problem still remains of finding the volume 
and some features are lost when using such calculations. 

References 

[1] Landau L. D., Lifshitz E. M., Fluid Mechanics, 1959, Pergamon Press, 
London. 

[2] Sedov L. L, 1959, Similarity and Dimansional Methods in Mechanics, 
Academic Press, New York. 

[3] Bisnovatyi-Kogan G. S., Silich S. A., Revs. Mod. Phys., 1995, V. 67, 
P. 661. 

[4] Klimishin I. A., 1984, Shock Waves in the Star Envelopes, Nauka, 
Moskow. 

[5] Hnatyk B. L, Pis'ma Astron.Zh., 1985, V. 11, P. 785. 

[6] Kompaneets A. S., 1960, Dokl. Akad. Nauk SSSR, V. 130, P. 1001 [Sov. 
Phys. Dokl, 5, 46]. 

[7] Silich S. A., Fomin P. I., Dokl. Akad. Nauk SSSR, 1983, V. 268, P. 861 
[Sov. Phys. Dokl., 28, 157]. 

[8] Korycansky D. G., Astrophys. J., 1992, V. 398, Pt. 1, P. 184. 



17 



[9] Kontorovich V. M., Pimenov S. F., Model of strong explosion in inho- 
mogeneous medium, Dopovidi NAN Ukrainy, 1996, No. 1, P. 54. (In 
English). 

[10] Kontorovich V. M., Pimenov S. F., Solar Physics, 1997, V. 172, P. 93. 

[11] Kontorovich V. M., Pimenov S. F., Model of strong explosion in a 
nonuniform corona — from solar radio bursts to relativistic jets, XXVIth 
Radio Astronomical Conference, Proc, St. -Petersburg, 1995, P. 197. 

[12] Kontorovich V. M., Pimenov S. F., II and IV types radio bursts in the 
strong explosion model in a nonuniform solar corona, Inform. Bulletin 
Astron. Ass. Ukraine, 1995, No. 7, P. 93. 

[13] Korn G. A. and Korn T. M., Mathematical Handbook for scientist and 
engineers, Mc Graw Book Company, N.Y., 1961. 

[14] Parker E. N., Astrophys. J., 1961, V. 133, P. 1014. 

[15] Hundhausen A. J. Coronal Expansion and Solar Wind, 1972, Springer, 
Heidelberg. 

[16] Sakurai A., Comm. Pure and Appl. Math., 1960. V. 13. P. 353. 

[17] Chernin A. D., Vistas in Astronomy, 1996, V. 40. P. 257. 

[18] Gvaramadze V. V., Pis'ma Astron. Zh., 1997. V. 23, P. 606. 

[19] Muhleman D. O., Anderson J. D., Astrophys. J., 1981, V. 247, P. 1093. 

[20] Lada C. J., Ann. Rev. Astron. Astrophys., 1985, V. 23, P. 267. 

[21] Kovalenko I. G., Shchekinov Yu. A., Astrofizika, 1985, V. 23, P. 363. 

[22] HST News. Supernova collision, http://www.stsci.edu (Usp. Fiz. Nauk, 
1997, V. 167, P. 778.) 



18 



